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Abstract 

Zero-mode waveguides (ZMWs) are confining light into attoliter volumes, enabling single molecule fluo¬ 
rescence experiments at physiological micromolar concentrations. Among the fluorescence spectroscopy 
techniques that can be enhanced by ZMWs, Forster resonance energy transfer (FRET) is one of the most 
widely used in life sciences. Combining zero-mode waveguides with FRET provides new opportunities 
to investigate biochemical structures or follow interaction dynamics at micromolar concentration with 
single molecule resolution. However, prior to any quantitative FRET analysis on biological samples, 
it is crucial to establish first the influence of the ZMW on the FRET process. Here, we quantify the 
FRET rates and efficiencies between individual donor-acceptor fluorophore pairs diffusing in aluminum 
zero-mode waveguides. Aluminum ZMWs are important structures thanks to their commercial avail¬ 
ability and the large literature describing their use for single molecule fluorescence spectroscopy. We 
also compare the results between ZMWs milled in gold and aluminum, and find that while gold has a 
stronger influence on the decay rates, the lower losses of aluminum in the green spectral region provide 
larger fluorescence brightness enhancement factors. For both aluminum and gold ZMWs, we observe 
that the FRET rate scales linearly with the isolated donor decay rate and the local density of optical 
states (LDOS). Detailed information about FRET in ZMWs unlocks their application as new devices 
for enhanced single molecule FRET at physiological concentrations. 

Introduction 

Forster resonance energy transfer (FRET) is one of the most widely used techniques in single molecule 
studies applied to life sciences [I]. FRET involves the nonradiative transfer of electronic excitation 
energy from an excited donor D* to a ground-state acceptor molecule A. The energy transfer efficiency 
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goes down with the inverse sixth power of the D-A distance, which makes FRET highly sensitive to the 
relative position of donor and acceptor fluorophores at the nanoscale. Therefore, FRET is often used 
to quantify the spatial relationship between two fluorophore-labeled sites in biological structures mm, 
study conformational changes in macromolecules m or follow interaction dynamics between proteins, 
DNA, RNA and peptide molecules W- 

FRET experiments are generally implemented on a diffraction-limited confocal microscope which 
requires pico to nanomolar concentrations so as to isolate a single molecule in the detection volume. 
To go beyond the restrictions imposed by diffraction, metal nanoapertures (also known as zero-mode 
waveguides ZMWs) are powerful tools for monitoring real-time dynamics of single molecules at micro¬ 
molar concentrations mm- By confining the illumination light to the bottom of the subwavelength 
aperture, detection volumes down to the attoliter range are commonly reached, reducing the detection 
volume by three orders of magnitude as compared to diffraction-limited confocal microscopy mm- 
Within ZMWs, experiments can be performed at the physiological concentration of molecules as found 
in living cells. ZMWs were successfully used on various investigations such as single-molecule DNA and 
tRNA sequencing P2HH], protein-protein interaction dynamics nans], and cell membrane receptors 
organization HZH2Q]. 

The use of ZMWs appears promising to extend single molecule FRET towards higher physiological 
concentrations. However, prior to any quantitative FRET analysis, it is crucial to establish first the 
influence of the ZMW on the FRET process. The question arises as FRET naturally competes with 
the donor radiative emission and the donor nonradiative energy losses to the metal structure, which are 
both known to be influenced by the photonic environment through the local density of optical states 
(LDOS) [2T|J22j . The LDOS is commonly defined as the number of electromagnetic modes per unit 
volume and frequency at the position of the dipole emitter where the energy can be released during 
the spontaneous emission process |22]. The energy release can occur by photon radiation or by local 
absorption and non-radiative transitions. The LDOS is therefore proportional to the emitter’s total 
decay rate (inverse of the fluorescence lifetime), including both radiative and non-radiative transitions. 

Several earlier works have considered the influence of photonic nanostructures on the FRET process: 
some conclude that the FRET rate depends linearly on the donor emission rate and the LDOS msM\, 
while some others report a FRET rate independent of the LDOS [281431] . Earlier experiments on 
aluminum C-shaped nanoapertures did not reveal noticeable changes of the FRET efficiency [32]. A 
similar trend was reported for circular aluminum ZMWs but without quantitative data m, while 
another study quantified a reduction of the FRET efficiency by 15% induced by the ZMWs j3lj- Recently, 
our group has investigated FRET in ZMWs milled in gold films, reported a linear dependence of the 
FRET rate on the LDOS and quantified a slight variation of the FRET efficiency with the aperture 
diameter [35] . 

Here, we build on our previous methodology [48] to investigate the role of aluminum ZMWs on the 
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FRET process between individual donor-acceptor fluorophore pairs on double stranded DNA linkers. 
ZMWs milled in aluminum are by far the most widely used structures as compared to other metals such 
as gold or copper [8 ti20ll32l - l34| . thanks to their commercial availability, low chemical reactivity and the 
availability of surface passivation protocols [36J. Due to the large interest raised by aluminum ZMWs 
and FRET separately, we believe that a clear quantification of FRET in aluminum ZMWs is necessary. 
Moreover, to provide a complete picture and discuss the role of plasmonic effects, we compare the results 
found with gold and aluminum, and find some significant differences on the fluorescence brightness and 
the fluorescence lifetime reduction notably. 

A key strength of our approach is the use of two different independent characterization methods 
to measure the FRET rate T fret and efficiency Efret- This is made possible by monitoring simul¬ 
taneously both the donor and the acceptor fluorescence photodynamics. The first method is based 
on the donor lifetime reduction in the presence of the acceptor. The second method is based on 
fluorescence burst intensity analysis accounting for the acceptor fluorescence increase in presence of 
the donor and the donor fluorescence quenching in presence of the acceptor. Both methods con¬ 
verge remarkably well towards similar results, as we show in Fig. 6 of the Supporting Information. 
The FRET efficiency is defined as the probability of energy transfer over all donor transition events: 
Efret = ^ fretda = ^ fret / (^ fret + ^ d) , where V D = T n,rad + ^ D,nonr is the total decay rate of 
the isolated donor (accounting for both radiative and non-radiative transitions) and T da = To + F fret 
is the total decay rate of the donor in the presence of the acceptor. As demonstrated for instance in [22], 
the LDOS is proportional to the isolated donor total decay rate T d- 

Results and Discussion 

In this work, we use circular nanoapertures with diameters from 110 to 400 nm milled by focused ion 
beam in a 150 nm thick aluminum or gold layer deposited on a glass microscope coverslip (Fig. [Hi). The 
FRET samples are provided by Atto550 donor and Atto647N acceptor covalently attached on double 
stranded DNA molecules with different distances between the donor and acceptor. To ease the discus¬ 
sion, we focus on two D-A separations of 20 base pairs (~ 6.8 nm) and 30 base pairs (~ 10.2 nm). The 
ZMWs are covered by a solution containing the FRET donor-acceptor pairs at around one micromolar 
concentration. The experiments monitor the fluorescence bursts generated as the fluorescent molecules 
cross the ZMW volume. Both the donor and the acceptor emission photodynamics are recorded simulta¬ 
neously with picosecond resolution in time-tagged time-resolved (TTTR) mode, resulting in fluorescence 
time traces as displayed in Fig. [T}x We have also checked that our conclusions are concentration invari¬ 
ant, see Fig. 8 of the Supporting Information. 

First we consider the effect of the ZMW on the fluorescence brightness of the isolated donor. Section 2 
of the Supporting Information details the use of fluorescence correlation spectroscopy (FCS) to quantify 
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Figure 1: (a) Scanning electron microscopy image of nanoapertures of 190 nm diameter milled in gold or 
aluminum with focused ion beam, (b) Typical fluorescence time traces for Atto550 donor and Atto647N 
acceptor diffusing in an aluminum nanoaperture of 190 nm diameter. The donor and acceptor molecules 
are fixed on a double stranded DNA linker with a donor-acceptor separation of 10.2 nm (30 base 
pairs). For each fluorescence burst exceeding the detection threshold in the donor or acceptor detection 
channel, a FRET efficiency is calculated (bottom trace). The binning time is 0.2 ms which is similar to 
the diffusion time of DNA samples in the ZMW (see Supporting Information Fig. 7). The total trace 
duration is 200 s. (c) Fluorescence brightness enhancement factors for the isolated donor for gold and 
aluminum apertures. Lines are guide to the eyes. 


the fluorescence brightness enhancement. Typical FCS correlation traces are shown in Fig. 7. We 
measure a typical fluorescence brightness enhancement of 3.2 in a 190 nm aluminum ZMW, while the 
fluorescence enhancement is only 2.35 in a gold ZMW of similar dimensions (Fig. QJ;). These results 
stand in good agreement with numerical simulations [37J38] and the behavior expected from previous 
measurements on red fluorescent dyes 11). al l. The difference between the fluorescence enhancement 
factors for aluminum and gold ZMWs stems mainly from the 550 nm laser excitation intensity which 
is lower in gold ZMW due to the increased losses for gold permittivity at 550 nm [37, 38], HU 01] (the 
phenomena leading to the observation of enhanced fluorescence are discussed with more details below). 
From the fluorescence intensity enhancement, it seems that aluminum ZMWs are therefore preferable 
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Figure 2: (a) Normalized donor fluorescence decay traces when no acceptor is present (circles) or when 
an acceptor is set at 10.2 nm distance (triangles) for the confocal reference and for a 190 nm aperture 
milled in aluminum or gold. Lines with symbols are numerical fits convoluted by the instrument response 
function (IRF). (b) Donor fluorescence lifetime as function of the aperture diameter for ZMWs milled 
in aluminum and for different donor-acceptor separations. Lines are guide to the eyes, (c) Same as (b) 
for apertures milled in gold. The reference lifetimes for the confocal setup are: 3.67 ns for the isolated 
donor, 3.13 ns for the donor with acceptor set at 10.2 nm, and 2.44 ns for the donor with acceptor set 
at 6.8 nm. All lifetimes are listed in the Supporting Information table S2. 


over gold for excitation wavelengths below 570 nm. Hereafter, we will investigate the influence of the 
metal on the fluorescence lifetimes and FRET rates in ZMWs. 

For each experiment, we construct the donor fluorescence decay trace by time-correlated single 
photon counting. Figure [2^i shows typical decay traces for the donor in confocal setup and in gold or 
aluminum ZMWs of 190 nm diameter. Two effects are readily seen on the raw decay traces. First, the 
donor emission dynamics are faster in the gold ZMW than in aluminum or confocal, as a consequence 
of the stronger influence of gold on the LDOS. Second, in all cases the donor emission dynamics are 
further accelerated by the acceptor presence, which opens a new decay channel by FRET. Consequently, 
the fluorescence lifetime of the donor in presence of the acceptor tda is shorter than the lifetime of the 
isolated donor td■ Numerical interpolation of the decay traces (see Experimental Section) quantifies 
the average donor fluorescence lifetimes for aluminum (Fig. Oo) and gold ZMWs (Fig. [5}:). A clear 
reduction of the donor emission lifetime (donor quenching) is observed as the D-A distance is reduced, 
with a similar trend observed for all nanoaperture diameters (Table S2 of the Supporting Information 
provides a list of all fluorescence lifetimes measured). We point out that an intrinsic limitation of our 
approach based on random diffusion is that it provides only the spatially averaged fluorescence lifetime, 
while lifetime variations are expected within the ZMW depending on the position and orientation of the 
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Figure 3: (a) Comparison of the decay rate enhancement (lifetime reduction) for the isolated donor as 
function of the ZMW diameter for gold and aluminum layers, (b) FRET rate enhancement as function 
of the ZMW diameter for different metals (aluminum: blue; gold: red) and donor-acceptor distance 
(6.8 nm: squares; 10.2 nm: triangles), (c) FRET rate T fret as function of the isolated donor decay 
rate F d for two donor-acceptor separations in ZMWs milled in aluminum (blue markers) or gold (red 
markers). 


emitter respective to the metal [431144] . Nevertheless, our net results remain relevant to describe the 
average properties of molecules randomly positioned inside ZMWs or randomly diffusing HU132H34]. 

The fluorescence lifetime td of the isolated donor decreases as the ZMW diameter is reduced, bringing 
the metallic walls closer to the dye (Fig. [2Jo,c) . Alternatively, the lifetime reduction can be seen as an 
increase of the isolated donor total decay rate T'p = 1 /to = T D,rad + r D,nonr, which accounts for both 
radiative and non-radiative transitions. Following the common definition in nanophotonics [22], is 
proportional to the LDOS at the donor emission wavelength. Figure [3ja, compares the enhancement 
of the isolated donor decay rate F d for aluminum and gold ZMWs. A clear difference is seen, with 
gold ZMWs providing ~ 2x faster photodynamics and higher LDOS, in agreement with numerical 
computations [51] and the presence of stronger plasnronic resonances for gold. 

In the presence of the acceptor, the total decay rate of the donor is increased by the additional FRET 
decay channel Tda = T d + r fret- The lifetime data in Fig. [2)o, c thus quantifies the FRET rate as 
r fret = r da ~Fd = 1/tda — 1/ td ■ For both aluminum and gold ZMWs, we observe similar increases 
of the FRET rate T fret and the isolated donor rate T ^ (Fig. [3[i,b), with a larger rate enhancement 
in the case of gold. Remarkably, our data indicates that the FRET rate T fret scales linearly with 
the isolated donor rate T^ for both metals and D-A distances (Fig. Eb). In the case of aluminum, the 
influence of the ZMW on the LDOS is less pronounced, and consequently the range of donor rates F d 
being probed is smaller than for gold. Nonetheless, the FRET rate in aluminum ZMWs follows a clear 
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Figure 4: Summary of the enhancement factors for gold and aluminum ZMWs. The total decay rate 
is for the isolated donor (no FRET). For both metals and D-A distances, the FRET rate enhancement 
evolves as the total decay rate of the isolated donor. The gains always refer to the confocal reference. 
See Experimental Section for details on the computation procedure. 


linear dependence with Tp, confirming the possibility to control FRET with nanophotonics [231148] . 

Another remarkable feature of FRET in ZMWs is the larger FRET rate enhancement in the case 
of increased D-A separations (Fig. [3)o) . For short D-A distances (on the order of the Forster radius 
or below), the direct dipole-dipole energy transfer dominates and the LDOS has a moderate effect on 
the FRET rate, especially for nanophotonic structures of the size of a half wavelength with limited 
field gradients m- However, for larger D-A distances and structures with more pronounced field 
confinement, a supplementary contribution from the energy transfer mediated by the nanostructure can 
further enhance the apparent FRET rate [231 [48]. 

To summarize the differences between aluminum and gold ZMWs, we quantify the different effects 
leading to enhanced photodynamics (Fig. [12] and Fig. 9 of the Supporting Information). As already 
seen in Fig. [Tfc, the fluorescence brightness appears significantly larger for aluminum than for gold. 
As a consequence of the increased losses for gold at the 550 nm illumination wavelength [31] [32]. the 
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Figure 5: (a) FRET efficiency histograms extracted from fluorescence burst analysis for 10.2 nm donor- 
acceptor separation in 190 nm nanoaperture milled in gold or aluminum. The FRET histogram found 
for the confocal reference is also displayed for comparison. Lines are Gaussian fits used to determine the 
center FRET efficiency, (b) FRET efficiency Efret as function of the aperture diameter for different 
metals (aluminum: blue markers; gold: red markers) and donor-acceptor distance (6.8 nm: squares; 
10.2 nm: triangles). The grey markers represent the values for the confocal reference. Lines are linear 
fits, (c) Distribution of FRET rates deduced from the efficiency histograms in (a) and the average 
isolated donor lifetime in Fig. [2b, c for 190 nm apertures and 10.2 nm donor-acceptor separation. 


excitation intensity enhancement is smaller for gold than for aluminum, which contributes to the smaller 
apparent fluorescence enhancement found for gold. Thanks to the stronger plasmon response of gold 
and probably also because of the oxide layer naturally present at the surface of aluminum, the radiative 
rate and total decay rate enhancements are ~ 2x larger for gold. For comparison, we add on Fig. [T2lthe 
values found for the FRET rate enhancement in the cases of 6.8 and 10.2 nm D-A distances. For both 
metals, the FRET rate enhancement is found in agreement with the isolated donor total decay rate. 
We also find that the FRET rate enhancement does not scale with the excitation intensity gain (as one 
would expect for FRET) nor with the total fluorescence brightness enhancement of the isolated donor. 
The case of the 190 nm aluminum ZMW is particularly illustrative: the fluorescence enhancement and 
local excitation intensity are high, but not the donor decay rate or the FRET rate. 

Complementary to the fluorescence lifetime analysis, our procedure also records the fluorescence 
bursts intensities from which the FRET efficiency Efret can be computed and collected in a histogram. 
The FRET efficiency is defined as the probability of energy transfer over all donor transition events: 
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Efret = ^ fret DA = ^fret/(Xfret + ^d) i where Tf, = T D,rad + ^D,nonr is the total decay rate of 
the isolated donor (without acceptor, accounting for both radiative and non-radiative transitions). The 
Experimental Section details the procedure to take into account the direct excitation of the acceptor 
by the laser light, donor emission crosstalk into the acceptor channel, and differences in the quantum 
yields and detection efficiencies of the donor and acceptor emission. Figure Or summarizes the FRET 
efficiency histograms for 10.2 nm FRET pairs in gold and aluminum ZMWs of 190 nm diameter and 
in the confocal reference case (non-normalized histograms are provided in the Supporting Information 
Fig. 9). Although the LDOS is clearly different between the three cases (Fig. [2^), all histograms have 
similar average values and widths. The trend is further evidenced by recording the average FRET 
efficiency as the ZMW diameter is varied (Fig. 0)). For the 6.8 nm D-A separation, we observe a 
small reduction of 6% of the FRET efficiency when the ZMW diameter goes down to 110 nm. The 
results for both aluminum and gold ZMW follow the same trend. For the 10.2 nm D-A separation, the 
FRET efficiency is almost constant for all ZMW diameters. The small 2% offset seen between gold and 
aluminum ZMWs in this case is an artefact related to the different sample solution in use, as confirmed 
by confocal measurements. Additionally, the FRET histograms can be used to compute the statistical 
distribution of the FRET rates Tfret = ^d Efret /(I — Efret ) using the separate measurements 
of the average isolated donor decay rates T^ (Fig. [2^i) and the FRET efficiency from burst analysis 
(Fig. Or). This leads to the distributions seen in Fig. [5fc, which illustrates the enhancement of the 
FRET rate T fret brought by the ZMW. 

Conclusions 

We have detailed the influence of aluminum ZMWs on the FRET process between individual donor- 
acceptor fluorophore pairs. We have measured separately the FRET rate and the FRET efficiency by 
measuring simultaneously the donor and the acceptor fluorescence photodynamics. We have found that 
both the FRET rate and FRET efficiency are consistent with the linear dependence of the FRET rate on 
the LDOS in ZMWs nanoapertures, and that accordingly, the FRET efficiency is marginally affected by 
the ZMW.We have also compared the results obtained with gold and aluminum ZMWs, and found some 
significant differences on the fluorescence brightness and the fluorescence lifetime reduction notably. 
While aluminum ZMWs have a lower influence on the decay rates than gold ZMWs of similar diameter, 
the lower losses of aluminum at 550 nm wavelength enable higher gains in local excitation intensity and 
fluorescence brightness. To compare between metals and find the most suitable ZMW material for a 
chosen illumination wavelength, the ratio —Re(e)/Im(e) of the real part over the imaginary part of the 
complex permittivity e is often used as a figure of merit [531 [55] . While gold has a high figure of merit 
in the red spectral region, aluminum performs better for illumination wavelengths below 560 nm [411142] 
(Fig. 13 of the Supporting Information). We anticipate that this quantification of the FRET process 
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in ZMWs will promote their application as new devices for enhanced single molecule FRET analysis at 
physiological micromolar concentrations . 

Experimental Section 

Nanoaperture ZMW fabrication. Nanoapertures are milled by focused ion beam (FEI Strata Dual 
Beam 235) on 150 nm thick aluminum or gold films deposited using thermal evaporation on standard 
150 fim thick microscope glass coverslips. 

DNA samples. Double-stranded DNA constructs of 51 base pairs length are designed with one 
Atto550 donor on the forward strand, and one Atto647N acceptor on the reverse strand. Fluorescently 
labeled and HPLC-purified DNA single strands are obtained from IB A (Gottingen, Germany), modified 
with the corresponding N-hydroxysuccinimidyl ester (NHS) donor and acceptor fluorophore derivatives 
of ATTO550 and ATT0647N. Fluorophores are covalently linked to an amino-C 6 -modified thymidine 
with NHS-chemistry via base labeling. The forward strand sequence is 

5’ CCTGAGCGTACTGCAGGATAGCCTATCGCGTGTCATATGCTGTTdCAGTGCG 3’. 

The reverse strand sequence is 

5’ CGCACTGAACAGCATATGACACGCGAT 20 AGGCTATCCT 30 GCAGTACGCTCAGG 3’. 

The distances between fluorescent labels are set so that the donor and acceptor are separated by 20 or 
30 base pairs (corresponding to ~ 6.8 and 10.2 nm separations respectively). The characteristic Forster 
radius computed for Atto550 and Atto647N in pure water is 6.5 nm. The reference sequences carrying 
only the isolated donor or acceptor are constructed with unlabeled complementary strand respectively. 
The strands are annealed at 10 fiM concentration in 20 mM Tris, 1 mM EDTA, 500 mM NaCl, 12 mM 
MgCl 2 buffer, and by heating to 95°C for 5 min followed by slow cooling to room temperature. Double 
stranded DNA stocks are diluted in a 10 mM Hepes-NaOH buffer, pH 7.5 (Sigma-Aldrich). 

Experimental setup. Experiments are performed on a confocal inverted microscope with a Zeiss 
C-Apochromat 63x 1.2NA water-immersion objective. The excitation is provided by a iChrome-TVIS 
laser (Toptica GmbH) delivering 3 ps pulses at 40 MHz repetition rate and 550 nm wavelength. We use 
on average 40 fiW excitation power to avoid saturating the fluorescent dyes. Filtering the laser excitation 
is performed by a set of two bandpass filters (Chroma ET525/70M and Semrock FF01-550/88). Dichroic 
mirrors (Chroma ZT594RDC and ZT633RDC) separate the donor and acceptor fluorescence from the 
reflected laser light. The detection is performed by two avalanche photodiodes (Micro Photon Devices 
MPD-5CTC with < 50 ps timing jitter) with 620 ± 20 nm (Chroma ET605/70M and ET632/60M) and 
670 ± 20 nm (Semrock FF01-676/37) fluorescence bandpass filters for the donor and acceptor channels 
respectively. The photodiode signal is recorded by a fast time-correlated single photon counting module 
(Hydraharp400, Picoquant GmbH) in time-tagged time-resolved (TTTR) mode. Each trace duration 
is typically of 200s. The temporal resolution for fluorescence lifetime measurements is 37 ps at half- 
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maximum of the instrument response function. 

Fluorescence lifetime analysis. The time correlated single photon counting (TCSPC) histograms 
are fitted using Levenberg-Marquard optimisation, implemented using the commercial software Sym- 
photime 64 (Picoquant GmbH). The model considers a single exponential decay reconvoluted by the 
instrument response function (IRF, as shown on Fig. [2ji). The time interval for fit is set to ensure that 
more than 85% of the detected count events are taken into account in the region of interest. As shown 
already in some of our earlier work [40)1451118 ] and by the results in Fig. [2ji, a single exponential fit is a 
satisfactory approximation to the fluorescence decay traces in nanoapertures (Fig. 11 of the Supporting 
Information). For gold aperture diameters below 250 nm, a background signal is detected due to the 
metal photoluminescence, which we take into account by adding a supplementary decay term with a 
fixed 5 ps characteristic time (see Fig. 11 of the Supporting Information). 

FRET efficiency analysis. For every detected fluorescence burst above the background noise, we 
record the number of photons in the acceptor channel n a and in the donor channel n d . For both the 
confocal and the ZMW configurations, the threshold for burst recognition is set to the sum of the mean 
plus one standard deviation of the summed trace of donor and acceptor channels. To avoid experimental 
artifacts in the FRET analysis, we also carefully characterize the optical response of the isolated donor 
and the isolated acceptor for each ZMW. In the case of the isolated donor, we measure the fraction a 
of photons from the donor that fall into the acceptor detection channel due to non-negligible spectral 
overlap between the donor emission and the acceptor detection window. For all our measurements, we 
find a constant a = 0.17 that is not affected by the ZMW. In the case of the isolated acceptor, we record 
the number nf® of photons that result from the direct excitation of the acceptor dye by the laser light. 
The FRET efficiency is then computed: 


Efret = 


n a - an d ~ n 


de 


n a - an d - n* + 7 n d 


( 1 ) 


Here 7 = r) a (j) a /r] d (j) d accounts for the differences in quantum yields (4> a and 4> d ) and fluorescence detection 
efficiencies (ry a and r] d ) between the acceptor and donor. For the confocal reference and the ZMWs, we 
estimate 7 = 1.3 in the case of our setup. The full trace analysis is implemented using the software 
Symphotime 64 (Picoquant GmbH). 

Quantification of the enhancement factors. The fluorescence enhancement r)F for the isolated 
donor is quantified by burst intensity analysis or alternatively FCS (Fig. [TJ^) [4UI[51j . Below the fluo¬ 
rescence saturation regime, three phenomena contribute to the fluorescence enhancement 77 ^: the gains 
in collection efficiency ry K , quantum yield ry^, and excitation intensity rj exc , so that r]p = ri K ri ( f ) r] exc . The 
local excitation intensity enhancement r] exc is deduced from numerical simulations by averaging the 
electric field intensity in the lower half part of the aperture [SUES]- The gain in collection efficiency is 
assumed to be equal to ?y K = 1.5 for both aluminum and gold 190 nm ZMWs, following earlier works 
on aperture directivity [451I5T] . The gain in quantum efficiency r] ( j > = r] rad /'i]t 0 t can be rewritten as the 
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ratio of the gain in radiative rate rj ra( i over the total (radiative + nonradiative) decay rate enhance¬ 
ment rjtot- The total decay rate enhancement rj to t corresponds to the lifetime reduction as measured 
in Fig. [3jr. The radiative rate enhancement can then be computed from all previous measurements as 
Vrad = VFVtot / {tlexcVK,) ■ The non-radiative rate (without FRET) enhancement is computed following the 
relation rjtot = 4>ilrad + (1 — < t>)Vnrad and the known quantum yield 0 = 0.8 for Atto550 in confocal case. 
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1. The two FRET measurement methods converge towards similar results 
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Figure 6: Comparison of the results obtained from analysis of the donor fluorescence lifetime (filled 
markers) and the acceptor/donor burst intensities (empty markers). Both methods converge remarkably 
towards similar results for both the FRET efficiency (a) and the FRET rate (b). The blue-grey solid 
regions indicate deviation of 8% from the average for D-A distances of 6.8 nm, and 10% for D-A 
distances of 10.2 nm. To extract the FRET rate from the burst intensity analysis, we use Tfret = 
T dE f ret /(I — Efret) where Tp = 1/td is the donor decay rate obtained from time-correlated lifetime 
measurements on the isolated donor case. The results shown here are for aluminum ZMWs. Similar 
results are obtained for gold nanoapertures [48j . 
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2. FCS quantifies the diffusion properties in ZMWs and the fluorescence brightness 
enhancement 


Fluorescence correlation spectroscopy (FCS) is used to quantify the number of molecules N contributing 
to the detected fluorescence signal and the mean translational diffusion time T d for the fluorophores to 
cross the detection volume [391150] . FCS performs a statistical analysis of the temporal fluctuations 
affecting the fluorescence intensity by computing the second order correlation of the fluorescence intensity 
time trace G^ 2 \t) = (5F(t).5F(t + T))/(F(t)) 2 , where SF(t) = F(t ) — (F) is the fluctuation of the time- 
dependent fluorescence signal F(t), r the delay (lag) time, and (.) stands for time averaging. The analysis 
of the FCS correlation function is based on a three dimensional Brownian diffusion model [491150] : 


G^ 2 \t) 


1 

N 



2 


1 + 


T 


l-T 


exp 



1 

(1 + rjT d )sJ\ + (1/k) 2 t/T d 


(2) 


where N is the total number of molecules, ( F ) the average intensity, ( B ) the background noise, T 
the fraction of dyes in the dark state, tj- the dark state blinking time, T d the mean diffusion time 
and k the aspect ratio of the axial to transversal dimensions of the analysis volume. The background 
noise (B) originates mainly from the back-reflected laser light and from metal photoluminescence. At 
40 gW excitation power, it is typically of 0.5 kHz, which is quite negligible as compared to the typical 
fluorescence brightness per molecule. Triplet blinking is quite weak for the Atto550 dye. The triplet 
fraction T converges to a value of 9% for the confocal reference, and vanishes for the FCS traces in 
ZMWs. Therefore, we decided to set T = 0 for the FCS analysis in ZMWs. The shape parameter k is 
fixed to k = 5 for the confocal microscope from the PSF calibration, and kept as a free parameter for 
the fits in the ZMWs. As already noted in our previous works [401151] . this parameter converges towards 
values close to 1 in ZMWs. 

To compute the FCS correlation, the fluorescence intensity trace is recorded on the Hydraharp400 
time-correlated single photon counting module (Picoquant GmbH) in time-tagged time-resolved (TTTR) 
mode. Each trace duration lasts typically 200s. The correlation is then computed using the software 
Symphotime 64 (Picoquant GmbH) with lag times ranging from 1 //s to 1 s. To avoid artifacts on the 
autocorrelation curve at sub-microsecond lag times related to photodiode afterpulsing, we implement 
the FLCS background correction as described in [52], The FCS curves are fitted using the model of 
Eq. (1) and Levenberg-Marquard optimization. This procedure quantifies the number of molecules 
N and their mean translational diffusion time r d . The brightness per molecule is then computed as 
Q = ((F) — ( B))/N . Comparing the brightness in the ZMW and the confocal reference then quantifies 
the fluorescence brighness enhancement factor rjp = QzMw/Qconfoc ■ Figure S2 presents typical FCS 
correlation data recorded in confocal diffraction-limited mode and in 120 nm and 190 nm aluminum 
nanoapertures. The fit results are summarized in Tab. SI. 
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Figure 7: Fluorescence intensity time trace and corresponding FCS correlation function for the confocal 
reference (a,b) and nanoapertures milled in aluminum (c,d). Thick lines in (b,d) are numerical fits 
using Eq. (1) model. The fit results are summarized in Tab. SI. The results for 190 and 120 nm 
ZMWs use the same concentration to enable a direct comparison between the fluorescence signals and 
detection volumes, (e) Normalized correlation functions showing the reduction of diffusion time across 
the detection volume as the ZMW diameter is reduced. 
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Confocal 

ZMW 190nm 

ZMW 120nm 

Number of molecules 

N 

1.2 

3.5 

1.3 

Diffusion time 

Td 

540 

170 

140 

Triplet fraction 

T 

0.09 

0 

0 

Blinking time 

t t (ms) 

40 

- 

- 

Shape ratio 

K 

5 

1.5 

1 

Average intensity 

(F) (kHz) 

6.8 

66.6 

13.2 

Background noise 

(B) (kHz) 

0.1 

0.5 

0.5 

Brightness per molecule 

Q (kHz) 

5.6 

18.9 

9.8 

Fluorescence enhancement 

VF 

- 

3.3 

1.7 

Concentration used 

C (ill) 

2 

1500 


Detection volume 

Veil (fL) 

1 

0.0038 


Volume reduction 

Ryol 

- 

265 



Table 1: The parameters obtained by fitting the FCS correlograms displayed in Fig. S2b,d are used to 
quantify the fluorescence brightness enhancement and the detection volume reduction. 
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3. The FRET results are invariant upon 4x concentration change 




N S 6 


N £ 2 


A/ = 0.4 


A/ = 0.2 


FRET efficiency [%] 


FRET efficiency |%| 


Figure 8: FRET histograms as function of the average number N of FRET pairs in the detection 
volume, as determined by FCS. The case IV ~ 1 is the one used for our experiments. No significant 
variation on the FRET histogram is seen when the concentration is raised up to 4-6 times. At very high 
concentrations (N ~ 16) deviations occur due to spatial averaging, yet this condition is always avoided 
in our study. The threshold for burst detection was set to the minimum value (mean of the trace) to 
better reveal small variations in the histograms. Consequently, the histogram widths are slightly larger 
here. 
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4. Comparison between FRET histograms 




Figure 9: Non-normalized experimental FRET histograms together with Gaussian fits used to determine 
the center FRET efficiency and the standard deviation. Events with apparent transfer efficiency below 
zero are also shown. We also display the histograms obtained on the DNA samples containing only the 
donor fluorophore to provide a reference for the zero FRET case. 
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5. FRET histogram widths 



O Confocal 
O Gold ZMW 
O Aluminum ZMW 
O Donor only 
-Shot noise limit ZMW 

♦ Deniz et al 

-Shot noise limit Deniz et al 


Figure 10: Standard deviation extracted from the Gaussian fits of the FRET histograms as function of 
the mean FRET efficiencies. For comparison to our data, we also report the values taken from Deniz et 
al Ref. |2]. Additional data on gold ZMW and confocal measurements are also taken from our previous 
work |48j . The shot noise limit is computed following the common simple approach o e = y/E( 1 — E)/T, 
as found in [21153] . E is the average FRET efficiency, T is the threshold level which is T = 20 for Ref. [2] 
and T = 22.2 for our data. 
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6. Fluorescence decay traces in gold have a supplementary fast decay related to gold 
photoluminescence 


a b 




Figure 11: Normalized donor fluorescence decay traces for a 190 nm aperture milled in gold (a) or alu¬ 
minum (b). Black lines represent the results of the exponential fit model reconvoluted by the instrument 
response function (IRF, grey trace). In the case of gold, a fast (~ 5 ps, limited by the resolution of our 
instrument) supplementary decay is visible, that we relate to the photoluminescence of gold excited at 
550 nm. 
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7. Table of donor fluorescence lifetimes measured in ZMWs 



Diameter (nm) 

Isolated donor 

td (ns) 

FRET 6.8 nm 

t da (ns) 

FRET 10.2 nm 

t da (ns) 

Confocal 

- 

3.67 

2.44 

3.13 

Aluminum 

120 

2.25 

1.52 

1.99 


190 

2.46 

1.72 

2.19 


295 

2.86 

1.96 

2.53 


385 

3.01 

2.05 

2.72 

Gold 


1.30 

0.90 

1.12 



1.65 

1.09 

1.42 


230 

1.97 

1.31 

1.67 


280 

2.08 

1.38 

1.83 


330 

2.26 

1.46 

1.92 


380 

3.02 

1.86 

2.50 


Table 2: Donor fluorescence lifetimes displayed in Fig. 2b,c. 
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8. Side by side comparison of enhancement factors for 190nm nanoapertures 


□ Donor fluorescence 

□ Excitation intensity 

□ Donor radiative rate 

□ Donor non-radiative rate 
■ Donor total decay rate 

□ FRET rate 10.2nm 


Gold 190nm Aluminum 190nm 



Figure 12: Summary of the enhancement factors for 190 nm ZMWs milled in aluminum or gold. The 
gain in collection efficiency is similar for both metals and amounts to 1.5x. 
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9. Choosing the most appropriate metal: a figure of merit as function of spectral 
range 



Figure 13: Figure of merit for plasmonic metals defined as the ratio of the absolute value of the real 
part over the imaginary part of the complex permittivity of gold (red) and aluminum (blue) [541155] . 
The metal permittivity data is extracted from [51. s 32|- The dashed line indicates the crossing between 
gold and aluminum near 560 nm. 
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